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ABSTRACT 

Markov  Chain  models  have  been  ased  to  forecast  stocks  in 
a  wide  range  cf  manpower  systems.  Studies  have  been  done  in 
many  areas  such  as  education  planning,  hospital  planning, 
manufacturing,  private  research  and  development,  a  women's 
military  unit,  the  civilian  work  force  supporting  the  U.S. 
Navy  and  a  state  police  organization.  This  study  looks  at 
such  systems  under  conditions  of  change  and  develops  the 
eguations  that  describe  the  steady  state  distribution  of 
personnel.  The  conditions  of  change  include  systems  where 
recruitment  is  constant,  increasing  (decreasing)  additively, 
or  increasing  (decreasing)  multi plica tively  and  systems 
where  the  changes  in  total  system  size  are  additively  or 
multiplicatively  increasing    (decreasing) .  Numerical    exam- 

ples utilizing  these  models  are  provided,  along  with  a 
computer  program  of  the  formulas  written  in  the  language 
APL. 
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I.  BACKGROUND 

A,   INTRODUCTION 

According  to  Bartholomew  and  Forbes  [Ref.  1:  p.  1]r 


Manpower  planning  is  often  defined  as  the  attempt  to 
match  the  supply  of  people  with  the  jobs  available  for 
them. 


This  is  a  concept  that  is  an  integral  part  of  military  force 
planning,  today.  All  branches  of  the  Department  of  Defense 
are  involved  with  studies  concerning  manpower  planning. 
Since  large  numbers  of  personnel  are  involved,  the  use  of 
manpower  models  to  determine  the  results  of  changes  to  the 
existing  system  before  they  are  actually  implemented  can 
prevent  costly  future  problems.  This  is  especially  true  in 
terms  of  having  the  right  number  of  personnel  in  the  rignt 
places  with  the  right  grade  levels  and  skills,  at  the  right 
time.  To  guote  again  from  Bartholomew  and  Forbes  [Ref.  1: 
P-  6"), 


A  manpower   model  is  a   mathematical  description   of  how 
change  takes  place  in  the  system. 


Much  work  has  been  done  in  the  area  of  manpower/ 
organizational  modeling.  The  text  by  Bartholomew  and  Forbes 
(referenced  above  and  hereafter  referred  to  as  BaF)  compiles 
information  en  the  subject  from  numerous  disciplines  and 
standardizes  notation  to  describe  some  of  the  basic  models 
and  statistical  technigues  used  in  manpower  planning.  A 
central  theme  of  this  work  is  the  idea  of  depicting  an 
organization  as  a  system  of  stocks  and  flows  [Ref-  1:  pp. 
3,4],   The  personnel  within  the  organization  under  study  are 


divided  into  classes,  or  categories,  by  specific  character- 
istics (for  example:  age,  time  in  service,  grade,  etc). 
The  numbers  of  personnel  within  categories  at  a  specific 
time   are     considered   stocks.  Flows   are      the    movements      of 

personnel  from  one  category  to  another  during  a  unit 
interval  of  time,  such  as  a  year- 
Various  mathematical  models  that  predict  future  stocks 
from  one  year  <or  other  time  period)  to  the  next  are 
described  in  EaF.  Of  particular  interest  here  is  the  tran- 
sition model  that  is  based  on  the  theory  of  Markov  Chains. 
In  addition  to  describing  the  model  [Bef.  1:  pp.  85-98],  BaF 
provides  a  BASIC  computer  program  [Eef.  1:  pp.  248-260],  to 
calculate  future  stock  vectors.  That  program  was  converted 
into  an  APL  program  at  the  Naval  Postgraduate  School  (NPS) 
and  is  available  on  the  NPS  mainframe  computer,  an  IBM 
system  370.  The  APL  program,  like  its  predecessor,  computes 
stocks  from  cne  year  to  the  next  based  on  input  such  as, 
initial  stocks,  recruitment  and  wastage  (or  attrition)  rates 
and  the  transition  probabilities  for  the  stocks.  This  can 
be  done  repetitively  and  theoretically  could  compute  future 
stocks    forever. 

There  are   two      versions    of   the    basic      Markov    Chain    model 
that    are     well    developed    in    BaF.  One   is    based    on      a    fixed 

flow  of  recruits  into  the  system  and  the  other  on  a  fixed 
total  system  size.  Also  of  interest  are  other  versions  (or 
options)  that  deal  with  a  system  undergoing  growth  where 
either  recruitment  or  total  system  size  grows  in  an  additive 
or      multiplicative    fashion.  The      APL      program    allows      the 

computation    of    future    stocks    for    the   cases    where: 

(1)  recruitment    is    fixed,      grows      additively    or    multi- 
plicatively; 

(2)  the    total      system    size    grows    additively      or    muiti- 
plicatively. 


Steady  state  is  the  condition  which  a  system  reaches, 
after  sufficient  time  has  passed,  so  that  the  stocks  no 
longer  change-  There  have  been  steady  state  models  devel- 
oped for  the  cases  where  there  is  a  fixed  flow  of  recruits 
into  the  system  and  where  the  total  system  size  is  fixed. 
In  both  these  cases  the  number  of  personnel  in  each  category 
remains  fixed  after  steady  state  is  reached  by  the  system. 
The  vector  whose  components  are  these  numbers  is  called  the 
steady  state  stock  vector  (S.S.S.V.). 

Under  conditions  of  growth,  the  steady  state  of  a  system 
must  be  redefined.  In  those  cases,  steady  state  does  not 
exist  in  the  above  sense,  because  the  system  is  undergoing 
constant  expansion  and  so  stocks  will  never  reach  a  fixed 
size.  However,  it  turns  out  that  the  stock  sizes  do  reach 
constant  proportions  and  so  it  is  possible  to  talk  about  a 
steady  state  distribution  vector  (S.S.D-V.).  Components  of 
this  vector  represent  the  proportions  of  personnel  that 
remain  in  each  category  forever  once  steady  state  has  been 
reached.  In  effect,  in  taose  cases  where  the  distribution 
of  personnel  among  the  categories  no  longer  changes  with 
time  the  system  has  reached  a  steady  state  "in 
distrib  ution. " 

Although,  at  the  present  time,  the  APL  program  at  the 
Naval  Postgraduate  School  predicts  future  stocks  for  a 
system  that  grows  under  conditions  mentioned  above,  it  is 
able  to  supply  a  steady  state  stock  vector  (S.S.S.V.)  only 
for  the  case  where  there  is  a  fixed  flow  of  recruits  into 
the  system-  The  only  way  to  get  the  steady  state  distribu- 
tion vector  under  the  various  conditions  of  growth  is  to 
continue  to  run  the  program,  i.e.,  to  forecast  stocks  for  as 
many  years  into  the  future  as  necessary  until  steady  state 
is  reached.  This  is  not  a  satisfactory  method  for  several 
reasons  : 


(1)  it  takes  too  long  to  run  the  model  when  the  system 
reaches  steady  state  only  after  many  years  have 
passed , 

(2)  it  is  a  trial  and  error  process  of  going  back  and 
forth,  in  time,  to  try  and  find  out  at  what  point 
in    time   steady  state    has   been   reached, 

{3)      occasionally,      even   if      the    program   is   run      out   as 
far      as      the      present     formatting      of      the      program 
allows,      there   can    still      be    instances    where  steady 
state   has   not   yet  been   reached- 
This  thesis      uses   the      * basic   prediction      equation* ,      as 
described   in    BaF  [Hef.    1:    pp.      86-88],    to  correct    the    afore- 
mentioned     deficiencies      by      developing      models      to      compute 
analytically      the      steady   state      distributions      for      systems 
under   varying   conditions    of    growth.  Also   developed    was   an 

API  program,  to  De  added  as  a  subroutine  to  the  main 
program,    that  computes  the  S.S.D.V.    for  these   options. 

The  next  section  in  this  chapter  sets  out  the  necessary 
notation  and  eguations  to  enable  the  reader  to  follow  the 
development  of  these  steady  state  models.  For  more  amplifi- 
cation on  the  casic  Markov  model  see  BaF  [Ref.  1:  pp  3-8, 
86-132]. 

B.       NOTATION     AND    BASIC    EQUATIONS 

Stocks,  representing  personnel  at  an  organization,  are 
divided    into    k   categories.       The   stock    vector 

2(t)    =     (    e1  (t)  ,r.2(t)  ,.    .     .    ,nK  (t)) 

is  a  row  vector  where  each  element  n  •  (t)  ,  represents  the 
number  of  personnel  in  category  i  at  time  t.  Time  is  meas- 
ured discretely  with  a  time  interval  of  unit  length  denoted 
(t-1,t) .  The  most  frequently  used  period  of  time  is  the 
year  although  it  is  sometimes  convenient,  or  necessary,  to 
use    quarters,    months,    etc. 
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Flows  describe      the    movements   of   personnel.  There   are 

two  basic  types  of  flows  in  a  system.  The  first  type  of 
flow  is  internal  to  the  system  and  can  be  thought  of  as 
transfers  from  one  category  to  another  within  the  system. 
Demotions  and  promotions  would  be  examples  of  this  type  of 
flow.  These  are  denoted  by  the  k  by  k  matrix,  N(t-1),  where 
each  element,  n(*:  (t-  1 ) ,  represents  the  number  of  personnel 
moving  from  category  i  to  category  j  during  the  time 
interval    (t-  1  ,t)  . 

The  other  type  of  flow  involves  the  two-way  transfer 
between  the  system  (organization)  and  the  outside  world. 
First,  there  are  wastage  (or  attrition  )  flows  of  the  number 
of  personnel  leaving  the  system  from  category  i,  during  the 
time  interval  (t-1,t),  denoted  by  n^K+.(t-1),  (i  =  1,2,-..,k 
and  t  =  1,2,...).  Then  there  are  recruitment  (or  accession) 
flows  of  personnel  into  category  j  from  the  outside  world, 
during  the  period  (t-1,t),  denoted  by  nD: (t) ,  (j  =  1,2,...,* 
and  t  =  1,2,...).  Following  BaF  [ Ref .  1:  pp.  3,4  ],  the 
notation  for  recruitment  flows  is  slightly  different  froir 
that      for    internal      flows   and      wastage    flows.  This   is      to 

remind  the  reader  that  in  most  organizations  attrition  and 
internal  flows  (transfers)  are  accounted  for  first  and  then 
recruitment    is    determined    at    the   end    of    the   time    period- 

The  Markov  Chain  model  uses  flow  rates  (or  protabili- 
ties)  rather  than  actual  numbers.  P,  a  k  ijy  k  matrix  of 
probabilities   has   elements 

lJ  nu    (t-1) 

There  are  certain  basic  assumptions  that  a  Markov  Chain 
model    needs    tc    meet    to   be   valid: 

(1)    the    transition   probabilities,       p..     ,      do   not  change 
over    time; 
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(2)    the    population     in    each   category  is      homogeneous   so 
that    p{«       represents   the   probability    of     each    indi- 
vidual  in    category    i      moving,      independently  of    any 
other  individual,   to  category    j. 
While  these      assumptions    may   not      always  be      met,      there 
have  been  enough      studies   done   using   the  model      to   show    they 
are   sufficiently   realistic  to    justify      the   use   of    the   model. 
Even   in   cases  where   these      assumptions    don't   hold    very    well, 
the   model   was   fcund   to  be    guite   useful   for    predicting   future 
stocks    [Refs.    2,3]. 

Next,  w  =  (w1  ,w2  ,...,wK)  ^s  defined  as  the  wastage 
(attrition)  rate  vector  where  each  element,  w^ ,  is  defined 
as 


w  •    = 

t 


n-  (t-1) 


Ncte    that, 


I 


p..      +    w      =    1       ,         i   =    1 ,2  ,..  .,  k 


■i-i     li 

Finally,      r   =    (r1  ,r2  ,.. .  ,  rK )       is   defined   as    the    recruit- 
ment   proportion    vector    where    each   element 

.        Ho,    (t) 


'J 


E  (t) 


is  the   proportion  of  total   recruitment  that  will   be  allo- 
cated to  category  j,  with 

K 

R<t)    =    2   no.    (t) 
j  =  i        J 
being      the      total      recruitment        during      the      time      interval 

(t-1,t).      This    implies      £  r-     =1. 

j=< 
BaF   [Ref.    1:    pp.       87,88],    snow    that   future    stocks   may    be 
calculated    using    the    eguation 
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n     (t)    =        >n;(t-1)p.,     ♦  E(t)r.     ,  (eqnl.1) 

fa*  I  J 

for  j=  1,2, ...rk  and  t  =  1,2,...  .  This  equation  states, 
that  the  stocks  in  category  j  at  time  t  are  equal  to  the 
number  of  personnel  that  moved  into  category  j  from  anywhere 
within  the  system  during  time  period  (t-1,t),  plus  the 
number  of  personnel  recruited  into  category  j  during  the 
same   period.      In   matrix    form    the   equation  is 

n(t)     =   n(t-1)P   +    R(t)r  t   =    1,2,...       .  (ecn    1.2) 

BaF   [fief.    1,:    p.    88]    refers    to   this  as    "the  basic    prediction 

eguation. " 

Another    eguation      dealing    with  changes  in      total   system 

size   is    also    developed    in    BaF   [fief.  1:    pp.  94-96].         To   see 

this,    note    that 

K 

N(t)     =     2nc(t) 
is    the    total    system   size    at    time    t    and 

«(t)     =   N(t)    -    N(t-1)  (eqn    1.3) 

represents  the  change  in  the  total  system  size  during  the 
time  interval  (t-1,t).  This  implies  that  total  recruitment 
during     (t-1,t)    can    be    expressed    as: 

K 

E  (t)     =      /n    (t-1)w.    ♦    M  (t)     >   0         .  (eqn    1.4) 

The  first  term  in  this  expression  is  recruitment  that  is 
done  to  replace  those  who  left  (attrition  from  the  various 
categories)  while  the  second  term  consists  of  recruitment 
that  is  due  to  the  chanqe  in  total  system  size. 

If  M(t)  =0,    then  only  those   personnel  wno   leave  the 
system  are   realaced,   which  implies   that  the   total  system 
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size  remains  fixed.  If  H  (t)  >  0,  new  personnel  are  supplied 
to  cover  an  increase  in  total  system  size,  i.e.,  to  fill  new 
jobs.  If  M  (t)  <  0,  then  some  jobs  must  be  eliminated  as 
soon  as  they  are  vacated.   Of  course,  E(t)  >  0  always. 

Referring   back  to   Eguation  1.1,    replace  R(t)    using 
Equation  1.  4  and  let 

q  ..  =  p.  +  w  .  r  •  . 
lJ  'J  L    J 

Then,  as  Ba F  [Ref.  1:  pp.  94,95],  shows 

K 

n-  (t)  =   yn-(t-1)q-.  +  M(t)r   ,  (egn  1.5) 

J  &  L         v  J 

for    j=    1,2, ,k    and    t   =    1,2,...       .       To   put    this    into   matrix 

notation  let 


Q   =    l%  }    . 


i, j    =    1  ,2, ,k    , 


or 


0_  =    P   +    w  'r    . 
Then    Equation    1.5    may    be    expressed    as 

n_(t)     =  n_(t-1)Q   +    M(t)r_      ,         t   =    1,2,...  {egn    1.6) 

Equations  1.2  and  1.6  form  the  basis  from  which  five 
steady  state  submodels  (or  options)  will  be  derived.  The 
next  section  introduces  these  options  and  the  direction  of 
development . 

C.   OBJECTIVES 

The  objective  of  this  thesis  is  tc  develop  steady  state 
equations  modeling  five  different,  conditions  of  growth. 
Each  of  these  conditions  will  be  addressed  separately  as  a 
"growth   option."     The   first   three   options    deal   with 
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different  types  of  recruitment.  They  use  Equation  1.2  where 
R  (t)  is  the  total  recruitment  during  the  time  period  (t-1,t) 
and  in  these  three  options,  R  (t)  is  the  only  term  changing. 
The  last  two  options  deal  with  changes  in  total  system  size. 
They  utilize  Equation  1.6  where.  M(t)  represents  the  change 
in  total  system  size  during  the  period  (t-1,t)  and  in  these 
two  options,  M(t)  is  the  only  term  that  changes.  A  complete 
list  of  these  five  options  is: 

(1)  Fixed  recruitment:   fi  (t)  =  Pf   t>1  ,   R>0   fixed, 
i.e.,  recruitment  is  always  constant; 

(2)  Additive  increase  (decrease)  of  recruitment: 
R(t)  =  R  +  (t-1) M  ,   t>1 

where  M  is   the  constant  amount  of   yearly  increase 
(decrease) ; 

(3)  Multiplicative  increase  (decrease)  of  recruitment: 
E(t)  =  Re*-i,  t>1, 

where  100(0-1)    is  the  percent  of   yearly  increase 
(decrease)  ; 

(4)  Additive  increase  (decrease)  in  system  size: 
M(t)  =  M   ,   t>1 

where  M  is   the  constant  amount  of   yearly  increase 
(decrease)  in  total  system  size; 

(5)  Multiplicative  increase  (decrease)  in  system  size: 

M(t)   =    (e-1)  e*-iN(0)        ,     t>i 

where    100(6-1)         is    the    percent   of      yearly    increase 

(decrease)  in  total  system  size. 
If  H  =  3  in  option  2  or  6  =  1  in  option  3,  then  those 
options  revert  to  option  1,  of  fixed  recruitment.  If  ?!  =  0 
in  option  4  cr  6  =  1  in  option  5,  then  M(t)  =  0  and  the 
result  is  a  system  with  a  fixed  total  system  size.  That 
case  will  be  dealt  with  as  subcases  of  options  <*  and  5-  The 
fixed  total  system  size  case  and  option  one,  where  R  is 
fixed,  are  the  only  cases  where  it  is  possible  to  acnieve  a 
steady    state      stock   vector,       i.e.,         where    the      stocks    don  * t 
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change  with  time.  In  all  other  cases  steady  state  will  be 
considered  to  have  been  reached  when  a  steady  state  distri- 
bution   of   personnel  is   achieved. 

Chapter  II  addresses,  in  depth,  the  analytical  develop- 
ment of  the  steady  state  equation  for  each  of  the  options. 
It  starts  with  an  account  of  the  work  done  on  option  1  in 
BaF-  Then  it  proceeds  with  the  derivation  of  the  steady 
state  equations  for  options  2  through  5  from  either  Equation 
1.2  or  Equation  1.6  .  It  also  presents  examples,  one  a 
four-grade  hierarchical  system,  the  other  a  three  grade 
non-hierarchical  system.  In  addition  to  formulating  analyt- 
ical models  fcr  options  2  through  5,  an  APL  program  has  also 
been  developed  that  will  compute  the  steady  state  distribu- 
tion lor  tnose  options.  This  program  is  contained  in 
Appendix  F.  Appendices  A  through  E  are  computer  printouts 
of  terminal  sessions  of  the  two  examples  for  each  of  options 
1    through    5. 
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II.  OPTION  DEVELOPMENT 

The  options  mentioned  in  Chapter  I  will  be  discussei 
here  under  the  usual  restrictions  of  manpower  systems. 
These  restrictions  imply  that  w- >  0  for  all  i=1,2,...,k, 
since  in  manpower  systems  attrition  must  be  allowed  from  any 
category-  This  means  that  the  p  matrix,  as  defined  in 
Chapter  I,  is  composed  of  non-negative  elements  and  has  row 
sums  all  strictly  less  than  one.  The  Q  matrix,  also  defined 
in  Chapter  I,  is  a  stochastic  matrix  (i.e.,  it  is  composed 
of  non- negative  elements  and  has  row  sums  all  equal  to  one). 
These  facts  along  with  the  following  two  mathematical  theo- 
rems are  needed  in  the  derivation  of  the  steady  state 
results  for  the  five  options. 

Theorem  _1 

If  P  is  a  matrix  composed  of  non-negative  elements  and 
row  sums  strictly  less  than  one  and  0  is  a  scalar  then 
the  matrix  6I-P  has  a  unique  inverse  for  all  values  of 
e    >   1    [Ref.    4:    p.    42    ]. 

Theorem    2 

If  Q  is  a  stochastic  matrix  and  6  is  a  scalar  then  the 
matrix  6I-Q  has  a  unique  inverse  for  all  values  of 
6    >    1    [  Ref.    4:    pp.    60,6  1    ]. 

A.       RECRUITMENT    OPTIONS 

The  following  three  options  are  all  variations  of 
Equation  1.2,  where  R  (t)  ,  representing  total  recruitment,  is 
the  term  that  takes  on  different  forms.  Consequently,  the 
steady    state    equations      for    these    options   are      also    obtained 
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in  somewhat  different  forms.     These  steady  state  equations 
are  derived  in  the  next  three  subsections. 

1 -   Option  J :   Fixed  Recruitment 

The  steady  state   equation  for  this  first  option  is 

developed  in   BaF  [Ref.  1:   p. 90].    The   analysis  presented 

here  is  a  brief  summary  of   that  in   BaF,   using  the  same 

assumptions.     Recalling   the  basic   prediction  equation, 

(Equation  1.2),  and  assuming  that 

limit   R  (t)  =  R    is  fixed 
t  — oo 

it  follows  that  the  steady  state  stock  vector, 

n_  =  limit   n  (t) 
t  —  oo 

also   exists    [ Bef .    4:       pp. 40-43,48-50  ].         Therefore,      letting 
t—  oo      in    Equation    1.2   results   in 

n  =    nP    +    Rr       , 

or 

n(l  -    P)    =  Rr      . 

Finally,   knowinq  the   matrix  I-p  has  a   unique  inverse  (see 
Theorem  1),  the  result 

n  =  Rr(I  -  P)-i  ,  (eqn  2.  1) 

is   obtained      for   tne    steady      state    stock      vector    (S.S.S.V.). 

After   steady    state    has   been    reached,    it    follows    that 

Mt)    =   n_     , 

for    "larqe"    t.        The   followinq   two      examples    will   be   used   to 
illustrate   the    use    of    this   steady   state   equation. 
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Example    (a) 

This  first  example  uses  data  from  Example  4.6  in  BaF 
[  Bef .  1:  p.  97],  which  is  taken  from  the  women  officers 
system  of  one  of  the  British  services.  It  is  a  hierarch- 
ical,   four    category  system. 

The   inital   stock  vector   is: 

n_(0)    =    (129,    74,    28,    11) 

The    transition   matrix    is: 

'.728  .10  2  0  0 

0  .83  .046  0 

0  0  .867  .033 

0  0  0  .902 


P   = 


The    total  recruitment   is: 

R(t)    =    35    ,       fixed    for   all    t>1 
The    recruitment    proportion   vector   is: 

r   =    (1,    0,    0,    0)       . 

Then 

.  272  -.  102  0  0 

0  .  17  -.046  0 

0  0  .133  -.033 

0  0  0  .098 


I   -   P    = 


and   the   inverse    is 


(I   -   P) 


—  i    = 


3.676  2.206  .763 

0  5.882  2.034 

0                0  7.519 

0                0  0 
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Using   Equation    2.1 


n  =    (35)  (1,    0,    0,    0) 


=     (128.66,     77.21,    26.705,    8.995)     . 

After   rounding   to   the   nearest   integer   the   S.S.S.V.    is: 

n   =     (129,    77,    27,    9) 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 
can  be  found  in  Appendix  A,  confirming  the  steady  state 
stock  vector  above.  As  the  printout  shows,  steady  state  is 
achieved  by  t  =  10. 

Example    (b) 

This   is    a    non-hierarchical,    three-grade   system. 

The    inital    stock   vector   is: 

n.{0)     =    (300,     200,    100) 

The    transition    matrix    is: 

60      .15      -05 
.75      .20        0 

05         0         .90 

The    total   recruitment    is: 

R(t)     =    100    ,     fixed    for    all    t> 1 
The    recruitment    proportion   vector    is: 

r   =    (.70,     .07,    .  23)     . 
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Then 


I  -   P    = 


and   the   inverse    is 

Z4.267         .8  2.133 

(I   -    P)-i    =  (      4  2  2 

\2.133         .4  11.066 

Dsing   Equation    2-1 

(4.  267         .8 
4  2 

2.133         .4 

=  (375.7,  79. 2,  417.8) 
After  rounding  to  the  nearest  integer  the  S.S.S.V.  is: 

n_  =  (376,  79,  418) 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 
can  be  found  in  Appendix  A ,  confirming  the  steady  state 
stock  vector  above-  As  the  printout  shows,  steady  state  is 
achieved  by    t    =    50. 

2-       Option    2:      Additive    Recruitment 

This    option      addresses   the      fixed   additive      increase 
(decrease!  in    recruitment      for    each      time    period.  Here, 

R   +     (t- 1)  fl,    replaces    R  (t)     in    Eguation    1.2    resulting    in 

£(t)     =n(t-1)P    +    Rr    ♦     (t-1)Mr_     .  (egn    2.2) 
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In    this  option,    M   >    0    must    be    true,      otherwise    there 
would   be  negative  recruitment    (R  (t)  <0)         which   does   not    make 
sense   in  manpower  models.      This   means    there  are   two   possible 
cases- 
Case    1:      M    =    0 

Then 

limit      B  (t)    =    E       , 

t—  oo 

which    is   identical    to    option    1. 
Case  2:         M    >   0 

Assume    that, 

limit     =     n_ 

t  — oo  tM 

exists-        This   is      the   same   as  assuming   that      for    "large"   t, 
n(t)  behaves    like      tMn_     .  Using    the      above      assumption. 

Equation   2.2    can   be   rewritten   as 

n  (t)  n(t-1)  B  (t-1)  M 

P       +  r      +     r 


tM                       tM  tM   ~~  tM 
or    - 

n(t)  n(t-1)  (t-1)                           B                    (t-1)M 

=  — m  •    p        +     r      +     r 

tM  (t-1)  H  t                              tM    _                   tM         — 

Taking    the    limit   as   t—  oo  results    in 

n  =      n  P  +      r 


or 


n"   (I-P)    =    r 


Since    the   I-P      matrix    has   a    unigue    inverse       (see    Theorem    1), 
the    result 

n*  =   r  (I-P)-i 
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is      obtained.  The      steady    state      distribution      vector     is 

achieved  by    normalizing  ft,    which   results   in 

n           £   [I-P]-*  (eqn    2.3) 

n    _  —  —  _  _ 

D1'  'BV 


After  steady  state  has  been   reached  it  follows  that 
the  formula 

n(t)  =  tMn    for  "large"  t, 

can  be  used  to  approximate  the  stocks  at  time  t.  The  two 
examples  from  section  A1  are  worked  out  below  to  illustrate 
Equation   2.3. 

Example    (a) 

This  is  a  hierarchical,  four  category  system. 

The  inital  stccK  vector  is: 

ji(Q)  =  (129,  74,  28,  11) 

The    transitioE   matrix    is: 

728       .102  0  0 

0  .83         .046         0 

P   = 

0  0  .867      .033 

0  0  0  .90  2 


The    total    recruitment    is: 

R  (t)    =  E   +     <t-1)M  where    E    =    35    and    M    =    5 

The    recruitment    proportion  vector   is: 

r   =    (1,    0,    0,    0) 
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Then 


I    -   P    = 


.  272 

-.  102 

0 

0 

0 

.  17 

-.046 

0 

0 

0 

.  133 

-.033 

0 

0 

0 

•     .098 

and   the   inverse    is 


(I   -   P) 


Using   Equation    2.3 


-l    = 


-"S 


„        I?  (1.    0,    0,    0) 

n  = 


-N 


ft1«  6.902 


=     (.53,    .32,    .11,    .04) 

is   the    steady   state   distribution   vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  100, 
900,  and  975  can  be  found  in  Appendix  E,  confirming  that  in 
steady  state  there  will  be  535?,  325?.,  11%,  and  4%  respec- 
tively, in  categories  1,2,3,  and  4.  As  the  printout  shows, 
steady  state  distribution  (to  two  decimal  accuracy)  is 
achieved  oy    t   =    975    . 

Example     (b) 

This   is   a   non-hierarchical,    three-grade   system. 

The    inital    stock   vector   is: 

n(0)     =    (300,     200,    100) 
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The   transition   matrix   is: 

60      .15      .05 
.75      .20        0 
05         0        .90 

The   total   recruitment   is: 

R  (t)    =   R    +     (t-1)M  where    R    =    100    and    M   =    25 

The    recruitment    proportion  vector   is: 

r   =    (.70,    .07,    .23)     . 

Then 

(.4  -.15      -.05 

-.75  .8  0 

-.05  0  .10 

and  the  inverse    is 

(4.267         .8  2.133 

4  2  2 

2.133         .4  11.066 

Using    Equation    2.3 


n*  (.7,    .07,    .23) 


L  =   ~ 


n1»  8.728 


=     (    .  43,     .09,     .43) 


is  the  steady  state  distribution  vector. 

The  computer   printout  of   the  terminal   session  for 
this  example,   showing   the  stocks  for  years   10,   50,   100, 
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350,  and  400  can  be  found  in  Appendix  B,  confirming  that  in 
steady  state  there  will  be  43$,  9%,  and  48%  respectively,  in 
categories  1,2,  and  3-  As  the  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t   =    400    . 


3-      Option   3:      Multiplicative   Recruitment 

This      option    addresses      the    multiplicative      increase 
(decrease)  in    recruitment      for   each      time   period-  Here, 

R8t_1    ,    replaces   Eft)     in    Eguation    1.2    resulting    in 

n(t)     =  n(t-1)P  +    E8t-ir        .  (egn    2.4) 

The  values  that  6  can  assume  play  an  important  part  in  the 
derivation  of  the  steady  state  equation.  First,  6  >0  must 
be  true,  otherwise  in  even  numbered  years  there  would  be 
negative  recruitment  (  B(t)<0  )  which  does  not  make  sense  in 
manpower  models-      There  are   three   cases   to    be   examined. 

Case    1     :         0<    8    <1 

Then 

limit      R  (t)    =    0       , 

t  — ►  oo 

because      et_l — ^0      when      0<    €   <1      .  Therefore,         letting 

t—  *>    in   Equation    1.2    results   in 

n    =   nP    , 

or 

n  (I   -    P)    =   0_-  . 

Since  I-P  has  a  unique  inverse  matrix  (see  Theorem  1) ,  the 
only  solution  to  the  above  eguation  is  n  =  p_.  This  means 
that   such   a    system    becomes  empty   in    steady    state. 
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Case  2:    6=1 

Then 


limit   E  (t)  =  R 

t—  oo 


which  is  identical  to  option  1. 

C  as  e  3  :    e  >  1 

Assume  that, 

i-  ■.   -(t)    - 
limit  —  =  n 

t  —  oo   e*   — 

exists.    This  amounts  to  assuming  that  for  "large"  t,   n(t) 
behaves  like   6*  r  .    Using  the  above  assumption,   Equation 
2.4  can  be  rewritten  as 

n  (t)        IL(t-1)  RG*-! 

•  P    + 


e*      e*-i  e        9*" 

Taking    the    limit   as   t  — *-  oo     results    in 


*       +>     1       „  R 

n  =    n«   -    «P      +        -    r 

~     9  9   "", 

or 

Ifr  [ei-p]    =  Rr     . 

Finally,      knowing   that    the   matrix   6I-P    has    a    unique    inverse, 
since   9   >    1     (see    Theorem    1)  ,     the    result 

d  =   Rr  [  9I-P  l~i 

is   obtained.         When    n    is    normalized,    the    resulting    vector   is 
n,    the    steady    state    distribution    vector   where, 


fr  Rr_[9I-P]-i  {eqn    2.5) 

n    = 


n1«  m 
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After  steady  state  has  been   reached  it  follows  that 
the  formula 

njt)  =  Q^i  for  "large"  t, 

can  be  used  to  approximate  the  stocks  at  time  t.    Two  exam- 
ples are  worked  out  below  to  illustrate  Equation  2.5 

Example  (a) 

This  is   a   hierarchical,    four  category  system. 

The    inital    stccK   vector   is; 

n(0)    =    (129,  74,    28,    11) 

The    transition    matrix   is.: 

728  .102           0           0 

0  .83         .0  46         0 

0  0            .867      .033 

0  0              0           .902 


P   = 


The    total  recruitment   is: 

fi(t)    =    Re*"-*       where    R   =    35    ,       9   =    1.0  1 
The    recruitment    proportion   vector    is: 

r  =    (1,    0,    0,    0)       . 

Then 


81    -    P   = 


and   the   inverse    is 
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[    ei  -  p]-i  = 


Using    Equation    2.5 


ft  (35) (1,    0,    0,    0) 


°        n1«  223.965 


=     (-554,     .314,    .101,    .031) 

is   the    steady  state   distribution    vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 
can  be  found  in  Appendix  C,  confirming  that  in  steady  state 
there  will  be  55%,  31%,  10f°,  and  3%  respectively,  in  catego- 
ries 1,2,3,  and  4.  As  tne  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t   =    50    . 


Example    (b) 

This   is    a   non-hierarchical,    three-grade   system. 

The    inital    stock   vector    is: 

£(0)     =    (300,     200,     100) 

The    transitioE    matrix    is: 

.60      -15       .05 
I   .75       .20         0 
.05  0         .90 
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The    total  recruitment   is: 

R  (t)    =   R©*-1    where   R   =    100    ,      ©    =    1-01 
The    recruitment    proportion   vector   is: 

r  =    (.70,    .07,    .23)     . 

Then 


ei  -  p  = 


.41 

-.  15 

-.05 

-.75 

.81 

0 

-.05 

0 

.  11 

and   the    inverse    is 


Using    Equation    2.5 


(4.03  -75         1.83 

3.73         1.92         1.69 
1.83  -34         9.92 


(.44,     .09,    .47) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 
can  be  found  in  Appendix  C,  confirming  that  in  steady  state 
there  will  be  44%,  9%,  and  47%  respectively,  in  categories 
1,2,  and  3.  As  the  printout  shows,  steady  state  distribu- 
tion (to  two  decimal  accuracy)  is  achieved  by  t  =  100  - 
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B.   TOTAL  SYSTEH  SIZE  OPTIONS 

The  options  in  this  section  are  variations  of  Equation 
1.6  where  HI  (t)  ,  representing  the  change  in  total  system 
size,  is  the  term  that  will  take  on  different  forms.  As  a 
result,  the  steady  state  equations  for  these  options  are 
also  obtained  in  somewhat  different  forms.  In  the  next  two 
subsections  the  steady  state  equations  will  be  derived  for 
the  system  size  options. 

1-   Option  4:   Additive  System  Size 

This  option  addresses  the  additive  increase 
(decrease)  in  total  system  size  for  each  time  period.  In 
this  case,  M  (t)  =  M,  a  fixed  amount  of  increase  (decrease) 
per  time  period-   Rewriting  Equation  1.6  results  in 

n(t)  =  n_(t-1)Q  +  Mr_  .  (eqn  2.6) 

There  are  3  cases  to  be  considered. 

Case  1:   M  <  0 

In  this  case  jobs  are  being  eliminated  by  a  fixed 
amount,  H,  each  year  which  implies  that  eventually  N(t)  =  0 
is  reached  in  a  finite  amount  of  time.  This  implies  that 
n  =  0_  is  the  steady  state  stock  vector.  This  means  that 
sucn  a  system  becomes  empty  in  steady  state. 

Case  2:   M  =  0 

When  M(t)  =  0,  N  (t)  is  fixed  and  can  be  denoted  N. 
It   follows    that    Eguation    1.6    can    be    written   as 

n  (t)     =   n_(t-1)Q      .  •  (eqn    2.7) 

Since   the 

limit        n(t)    =    n 
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exists  [Ref.  4:  pp.  40-43,  48-50  ],  Equation  2.7  can  be 
written   as 

n   =    n   Q      , 
or 

n_    (I    -   Q)    =   0      . 

Because  the  first  column  of  I-Q  is  a  linear  combination  of 
all  the  other  columns  the  inverse  of  I-Q  does  not  exist. 
So,  to  solve  for  n,  it  is  necessary  to  introduce  the  addi- 
tional constraint  of  n1_*  =  N  [Ref.  5:  p.  80  ].  Replacing 
the  first  column  of  the  I-Q  matrix  with  a  column  of  ones, 
denoted  by  (I-Q) *,  and  replacing  the  first  element  of  £  with 
N,    denoted    0>  *  ,    results    in    the    following   equation 

_n     (I-Q)  *   =    0* 

Taking    the    inverse    results  in 

_n    =    0_*[  (I-Q)*]-*         ,  (egn    2.8) 

the      steady    state      stock      vector.  After    steady      state      is 

reached,  then  n_(t)  =  n.  for  "large"  t.  Two  examples  are 
worked    cut    below    to  illustrate   Equation    2.8    . 

Example    (a) 

This   is    a   hierarchical,    four    category    system. 

The    inital    stccK  vector   is: 

n(0)     =     (129,    74,    28,    1  1) 

The    transition    matrix    is: 


P   = 


728       .102  -     0  0 

0           .83  .046  0 

0              0  .867  .033 

0              0  0  .902 
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r 


and 

(.898  .102  0  0 

.124  .83  .046         0 

.1  0  .867  .033 

.098         0  0  .902 

The    0*    vector   is: 

0*    =     (    242,    0,    0,    0    )     . 
The    recruitment    proportion  vector   is: 

£»(1,    0,    0,    0)        . 

Then 


I  -   Q    = 


Replacing  the  first  column  with  ones  gives 


(I  -  Q) 


and    the   inverse    is 


*    = 


1 

-.  102 

0 

0 

1 

.  17 

-.046 

0 

1 

0 

.  133 

-.033 

1 

0 

0 

.  098 

((I    -    Q)*) 


—  l    = 


.533 

.320 

.111 

.037 

-4.582 

3.133 

1.084 

.  365 

-5.354 

-3.212 

-    6.408 

2.158 

-5.435 

-3.261 

-1.128 

9.  824 

Using    Equation    2. 
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.533  .320  .111 

4.582  3.133  1.084 
n   =     (24  2,    0,     0,     0)  , 

-5.354  -3.212  6.408 

-5.435  -3.261  -1.128 


=     (128.99,     77.44,    26.86,    8.95)     . 

After  rounding  to  the   nearest   integer   the   S^S.S.V.    is: 

n_  =     (129,    77,    27,    9)     . 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100, 
can  be  found  in  Appendix  D,  confirming  the  steady  state 
stock  vector  above.  As  the  printout  shows,  steady  state  is 
achieved  by  t  =  50 

Example    (b) 

This   is   a    non-hierarchical,"  three-grade   system. 

The    inital    stcck   vector   is; 

n(0)    =    (300,     200,    100) 

The    transition    matrix    is: 

(.60       .15      .05 
.75       .20         0 
.05         0         .90 

and 

(.74  .164  .096 

.785  .2035  .0115 

.085  -0035  .9115 

The    0*    vector  is: 

0*   =    (   600,    0,    0    )     . 
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The    recruitment    proportion  vector   is: 

r   =    (.70r     .07,    .23)     . 
Then 

(.26  -.164  -.096 
-. 785  .797  .  -.012 
-.085      -.004         .089 

Replacing  the  first  column  with  ones  gives 

(1  -.164  -.096 
1  .797  -.012 
1    -.004    .089 

and   the    inverse    is 

(.430  .090         .479 

-.611       1.127      -.516 
-4.888       -.981       5.8b9 

Using   Equation    2.8 

(.430  .090  .479 

-.611  1.127         -.516 

-4.868         -.981         5.869 

=     (    258.0,     54.4,    287.4)     . 

After   rounding    to   the    nearest    integer    the    S.S.S.V.    is 

n   =     (258,    54,    237)     . 
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The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50  and  100r 
can  be  found  in  Appendix  D,  confirming  tne  steady  state 
stock  vector  anove.  As  the  printout  shows,  steady  state  is 
achieved  by  t  =  50  . 

Case  3:   M  >  0 

Equation  1-6  can  then  be  rewritten  as 

n(t)    n  (t- 1)         M 

Q   + 


N  (t)      N(t)  N  (t)  -  . 

Iterating  the  eguation,   N(t)  =  N(t-1)  +  M,   for  t  =  1,2, 

will  shew  that,   N  (t)  =  N  (0)  +  tM    for  all  t  =  1,2,...   and 


therefore, 


M  N(t-1) 

limit    =0  and        limit    =    1 

t— ••      N(t)  t  —  <*        N(t-1)     +    B 


Then   writing 


n_(t)  Jl(t-1)    N(t-1) 

=    • 


N  (t)  N  (t-1)     [N(t-1)     +    H]  N(t) 

and    taking    the    limit    as   t — +■  oo     yields 

n    =   J  Q    ,    or      n_    (I-Q)    =    0_ 
where 

n(t) 


n    =    limit 


t— oo       N  (t) 

is   the    steady    state    distribution    vector. 

This  is  identical  to  the  M  =  0  case,  except  here  the 
constraint  n1_*  =  1  must  hold.  Therefore,  by  replacing  the 
first      column   of      the      I-Q  matrix      with      a      column    of      ones. 
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denoted    by     (I-Q)*,    and  replacing    the   first    element   of   0    with 
1r    denoted    0_**,    the   solution    is 

n.  =    0_**[  (I-Q)  *]-!       ,  (eqn    2.9) 

the    steady   state  distribution    vector- 
After  steady  state  has  been      achieved,      _n(t)       can    be 
found   using 

njt)     =    Q[N(0)     +    tfl]n 

for  "large"  t-    Two  examples  are  worked  out  below  to  illus- 
trate Eguation  2.9  . 

Example  (a) 

This  is  a  hierarchical,  four  category  system. 

The  inital  stock  vector  is: 

n.(0)  =  (129,  74,  28,  11) 

The    transition   matrix   is: 

'.728  .102  0 

0  .33  .046 

0  0  .867 

0  0  0  .902 


P   = 


and 


Q   = 


.898  .102  0  0 

.124  .33  .046         0 

.  1  0  .867  .033 

.098  0  0  .902 


The    0**    vector    is: 

0**    =     (    1,    0,    0,    0    ) 
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The    recruitment    proportion  vector   is: 

r  =    (1,    0,   0,    0)       . 
Then 


I   -    Q   = 


Replacing   the  first   column  with   ones,    gives 


(I    -   Q)*    = 


and   the  inverse    is 

.533  .320  .111 

-4.582  3.133  1.084 

-5.354  -3.212  6.408 

-5.435  -3.261  -1.128 


(d   -    Q)*) 


—  i    = 


.533 

.320 

.111 

.037 

-4.582 

3.133 

1.084 

.365 

-5.354 

-3.212 

6.408 

2.  158 

-5-435 

-3-261 

-1.128 

9.  824 

Using    Eguation    2.9 


n  =    (1,    0,    0,    0) 


=     (.53,    .32,    .11,    .04) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  100, 
900,  and  950  can  be  found  in  Appendix  D,   confirming  that  in 
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steady  state  there  will  be  53%,  32%,  11%f  and  4%  respec- 
tively, in  categories  1,2,3,  and  4.  As  the  printout  shows, 
steady  state  distribution  (to  two  decimal  accuracy)  is 
achieved  by  t  =  950 

Example  (b) 

This  is   a   non-hierarchical,    three-grade  system. 

The    inital   stcck  vector  is: 

n_(0)    =    (300,     200,    100) 

The    transition   matrix    is: 

(.60       .15       .05 
.75       .20         0 
.05         G         .90 

and 

(.74  .164         .096 

.785         .2035       .0115 
.085         .0035       .9115 

The    0**    vector    is: 

0**    =     (    1,    0,    0,    )     . 

The   recruitment    proportion  vector   is: 

r   =    (.70,    -07,    .23)     . 
Then 

(. 26  -. 164  -.096 
-.785  .797  -.012 
-.085      -.004         .089 

Replacing   the   rirst    column   with   ones,    gives 

-.164      -.096 

.797      -.012 

-.004         .089 
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and   the   inverse    is 


.430         .091         .479 
-.611       1.127      -.516 
k-4.888      -.98*1      5.869 


-430 

.091 

.479 

-.611 

1.  127 

-.516 

-4.888 

-.981 

5.869 

Using    Equation    2.9 


n  =     (1,0,     0   ) 


=  (  .  43.  .09,  .48) 
is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  100 , 
350,  and  400  can  be  found  in  Appendix  D,  confirming  that  in 
steady  state  there  will  be  435?,  9£,  and  48%  respectively,  in 
categories  1,2,  and  3.  As  the  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t  =  400   . 

2-   Option  5:   Multiplicative  System  Size 

This  section  addresses  the  option  of  multiplicative 
increase  (decrease)  in  total  system  size  for  eacn  time 
period.   Starting  with 

N  (t)  =  eN(t-1),   for  t  =  1,2, 

it    follows    that 

N  (t)    =    e*N(0)#       for   t    =    1,2,.-. 
from    which    M  (t)     can    be   determined   as 

M(t)    =    N(t)     -    N(t-1)       =        (6-1)  e^-iN  (0)     . 


40 


Then,       letting    (6- 1)  9r_lN  (0)       replace    M(t)       in    Equation     1.6, 
results   in 

n(t)     s   n(t-1)Q  ♦    [(G-DS^JNIO)]   r      -  (eqn    2.10) 

First,      6   >    0   must    be    true  otherwise    negative    or   zero   system 
size    may  result.      This    leaves    two   cases   to    be    examined. 

Case    1:  0    <   6  <    1 

In    this   case 

limit   N  (t)     =0 

t  — »■  <so 

because   as    t — *-  oo     ,    et- l — *■  0  for    0    <    6   <1    .         Finally,      if 
©   =    1    then    9-1    =   0    and    N  (t)    =0    for    all   t.       Therefore, 

n_  =   limit   n(t)     =  0 

t  -»  oo 

must    also    hold.      This    means    that    such    a   system    becomes    empty 
in   steady   state. 

Case   2:  9    >    1 

Equation    2.10    can    be    rewritten    as 

n(t)  n(t-1)  (9-1)9"t'-iN(0) 

q      +     r 


N  (t)  N(t)  N(t) 

Since   N  (t)     can    fce   written    as    9N(t-1)    or   as    9**N(0)  ,    the    above 


equation   gives 


n_(t)  n.(t-1)       1  (9-1)  e*-iN(0) 


-    .q      +  — - — — r 


N  (t)  N(t-1)       9  brN  (0) 

or 

n  (t)  n  (t-1)       1  9-1 

=    •  -   «Q      *     r 

N  (t)  N  (t-1  )        9  9      ~~ 

Takinq    the    limit   as    t — »    oo    results    in 
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~  1  <e-i) 

n      =      n   -   Q      +    r 

-  -  G  9        — 


or 


n   [GI-Q]      =      (6-1)r 

where 

n  (t) 

n    =    limit     

—        t-+  «o  N(t) 

is  the    steady  state   distribution    vector- 
Since   6    >    1      in   this      case,       the      matrix    6I-Q      has    a 
unigue    inverse   according    to    Theorem    2,    so   the   result 

n    =    (0-1)  r[  ©I-Q]-*  (egn    2.11) 

is   obtained    for    the  steady   state    distribution    vector.       After 
steady    state    has   been    reached    it    follows   that    the    formula 

n(t)    =   8trN  (0)     n 

can    be    used    tc    find    the   stocks   at    time    t.         Two    examples    are 
worked    out    below   to   illustrate   Equation   2-11    . 

Example    <a) 

This   is   a   hierarchical,    four    category    system. 

The    inital   stock   vector   is: 

n<0)     =    (129,    74,    28,    11) 

The    transition    matrix    is: 

.728       .102         0  0 

0  .83      ".046         0 

P   = 

0  0  .867      .033 

0  0  0  .902 
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and 


.  898 

.102 

0 

0 

.  124 

-83 

.0  46 

0 

.  1 

0 

.867 

.033 

.098 

0 

0 

.902 

Q  = 


The    multiplicative   factor   is: 

9      =    1.01       . 

The    recruitment    proportion   vector  is; 

r  =    (1,    0,    0,    0)       . 

Then 


.112  -.  102  0  0 

-. 124  . 18  -.046  0 
GI    -    Q    = 

-.  1  0  .143  -.033 

-.098  0  0  .108 


and    the   inverse    is 


(91  -    Q) 


—  i    = 


55.4  1  "31.40  10.  10  3.09 

51.04  34.48  11.09  3.39 

50.35  28.53  16.  17  4.94 

50.  28  28.49  9.  17  12.06 


Using    Equation    2.11 


n   =     (.01)     (1,0,0,0) 


55.41  31.40  10-10 

51.04  34.48  11.09 

50.35  28.53  16.  17 

50.28  28.49 


=     (-554,     .314,    .101,    .031) 
is   the    steady   state   distribution    vector. 
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The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  and  100, 
can  be  found  in  Appendix  E,  confirming  that  in  steady  state 
there  will  be  55^,  31%,  10%,  and  3%  respectively,  in  catego- 
ries 1,2,3,  and  4.  As  the  printout  shows,  steady  state 
distribution  (to  two  decimal  accuracy)  is  achieved  by 
t  =  50   . 

Example  (b) 

This  is  a    non-hierarchical,    three-grade  system. 

The    inital    stcck    vector   is: 

n(0)    =    (300,     200,    100) 

The    transition   matrix    is; 

/.60      .15      .05 
P    =       [    .75      .20        0 
\  .05         0         .90 

and 

(.74     .164    .096 
.785    .2035   .0115 
.085    .0035   .9115 

The  multiplicative  factor  is: 

e  =  1.01     . 

The    recruitment    proportion   vector   is: 

r   =    (.70,    .07,    .23)     . 

Then 


ei    -    Q    = 


.27        -.164         -.096 
-.785         .8065      -.0115 
-.085      -.0035         .3985 
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and   the   inverse    is 


4  5.32 

9.41 

45.  27 

44.69 

10.52 

44.79 

40.70 

8.50 

50.  81 

45.32 

9.41 

45-27 

44.69 

10.52 

44.79 

40.70 

8.50 

50.81 

(91   -    Q)-»    = 

Osing    Equation    2.11 
n  =     (.0  1)     (.7, .07,  .23) 

=     (.44r    -09r    .47) 

is  the  steady  state  distribution  vector. 

The  computer  printout  of  the  terminal  session  for 
this  example,  showing  the  stocks  for  years  10,  50,  and  100, 
can  be  found  in  Appendix  E,  confirming  that  in  steady  state 
there  will  oe  44??,  9%,  and  47^  respectively,  in  categories 
1,2,  and  3.  As  the  printout  shows,  steady  state  distribu- 
tion (tc  two  decimal  accuracy)  is  achieved  by  t  -  100  . 
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Ill-  SOMIAEI 

This  chapter  is  a  synopsis  of  the  options  that  describe 
the  five  different  conditions  of  growth.  In  all  options, 
the  restriction  of  wt  >  0  for  all  i  =  1,2,...,k  was  assumed, 
since  attrition  must  be  allowed  from  any  category  in  a 
manpower  system.  Options  1  through  3  employed  the  basic 
prediction  eguation  (Eguation  1.2)  and  used  constant,  addi- 
tively  and  multiplicatively  increasing  (decreasing)  recruit- 
ment, R(t),  to  derive  the  steady  state  behavior  of  a 
manpower  system.  For  options  4  and  5,  Eguation  1.6  used 
additively  and  multiplicatively  increasing  (decreasing) 
total  system  size,  M(t),  to  derive  the  steady  state  behavior 
of  a  manpower  system. 

Derivation  of  the  steady  state  stock  vector  or  distribu- 
tion vector  required  some  standard  results  of  matrix  algebra 
summarized  in  Theorems  1  and  2  of  Chapter  II.  In  addition 
to  the  analytical  derivation  of  the  steady  state  eguations, 
examples  have  been  worked  out  in  each  case  (excepting  some 
trivial  cases)  to  illustrate  the  use  of  the  analytical 
results.  The  formulas  have  been  programmed  in  the  AFL 
language  and  the  program  listing  is  given  in  Appeniix  F. 
This  program  (called  a  function  in  APL)  has  been  integrated 
into  the  main  APL  program  at  NPS  which  forecasts  future 
stocks  of  a  manpower  system  using  the  Markov  Chain  theory 
summarized  in  Chapter  I.  Appendices  A  through  E  show  the 
printouts  of  the  terminal  sessions  where  the  examples 
referred  to  above  are  worked  out  to  verify  the  steady  state 
results  obtained  in  each  analytical  example. 
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APP ENDIX  A 
COMPUTER  PRINTOUTS  OF  EXAMPLES  FOR  OPTION  1 

1.   Example  (a) 

START 

DO  YOG  WISH  TC  ENTER  DATA? 

0  NC 

1  YES 

1 

ENTER  THE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV    HIERARCHICAL 

2  MARKOV    LENGTH    OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 

1 
ENTER    N  (INITIAL    STOCK    VECTOR) 

129    74    28    11 

ENTER    THE    PROMOTION    RATE    VECTOR. 

THIS    VECTOR    SHOULD    INCLUDE   THE   PROMOTION    BATES    FOR    THE    FIRST 
3    CLASSES. 

.  102    .046    .033 
ENTER    THE    WASTAGE    RATE    VECTOR. 

THE    4TH    VALUE    SHCULD    INCLUDE    THE    TOTAL    RATE    FOR       THAI 
STATE    DUE    TO     EITHER    PROMOTION    OR    WASTAGE. 

.  17    .  124    . 1    .098 

ENTER  THE    NUMEER   OF    THE    RECRUIT    TYPE 

1  FIXED    RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM     SIZE) 

1 
ENTER    R     (RECRUITMENT    VECTOR) 

35   0    0    0 

ENTER    THE    PERCENT    CODE 

0  NO    GRADE    PERCENTAGES 

1  GRADE    SIZE    AS    PERCENT    OF    TOTAL    GRADE    SIZE 

2  GRADE     SIZE    AS    PERCENT    OF    ORIGINAL    GRADE    SIZE 
7       QUIT    PROGRAM 

1 

WOULD    YOU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 
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P    MATRIX 
0.728    0.102    0  0 

0  0.83       0.046    0 

0  0  0.867    0.033 

0  0  0  0.902 

N    VECTOR 
129    74    28    11 
OPTION    =1 

R    VECTOR 
35    0    0    0 

WOULD    YOU    LIKE    TO    CHANGE    ANY    OF    THE    DATA? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 
ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
1  0 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 
T 

N 

PERCENT 

R 

0            1 
2 
3 

4 
TOTAL 

129     I 
74     i 
28    | 
11     | 

242    { 

t           "| 

!           100) 

10            1 
2 

3 
4 
TOTAL 

129    | 
77     i 

27    | 
10     | 

242    I 

53) 

32) 

"I 

t                 100) 

35 


DO    YOU    WISH    TO    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PRCGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
50 

DO    YOU     WISH     TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 

50            1  129  (                     53) 

2  77  32) 

3  27  (  11) 

4  9     {  4)  - 

TOTAL  242     (  100)  35 
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DO    YOU     WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

1 

ENTER    THE    NUMBER    OF    THE    YEAR     YOU    WISH    TO    SEE 
100 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 

100     1 

129  | 

'        53) 
32 

2 

77  I 

3 

27  I 

11 

4 

9  I 

V 

TOTAL 

242  | 

*                 100 

35 


DO    YOU    WISH    TO    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 


DO    YOU    WISH    TO    SEE    THE    STEADY    STATE    VECTOR? 

0  NO 

1  YES 


1 

999     1 
2 

3 
4 
TOTAL 

129  i 
77  i 
27  | 

9  I 

242  | 

[        53) 

32 
11 

100) 

ARE  YOU  THROUGH? 

0  NC 

1  YES 

35 


Examj^Ie     (b) 


START 

DC    YOU    WISH    TC    ENIER    DATA? 

0  NC 

1  YES 

1 

ENTER  THE    NUMEER    OF    THE    MODEL    TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    LENGTH    OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 
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ENTER    N  (INITIAL    STOCK    VECTOR) 

300    200    100 

ENTER    P     (TRANSITION    MATRIX)     BY    ROWS 
ENTER    1  TH    ROW 

.6.15    .05 

ENTER    2IH    ROW 

.75    .2    0 

ENTER    31H    ROW 

.  05    0    .9 

ENTER  THE    NUMBER   OF    THE   RECRUIT    TYPE 

1  FIXED    RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 
H  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 

1 

ENTER    R     (RECRUITMENT    VECTOR) 

70   7    23 

ENTER    THE    PERCENT    CODE 

0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PEP.CENT  OF  ORIGINAL  GRADE  SIZE 
7   QUIT  PROGRAM 

1 

WOULD    YOU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 

P  MATRIX 
0.6   0.15  0.05 
0-  75  0.2   0 
0.  05  0     0.S 
N  VECTOR 
300  200  100 
OPTION  =1 

R  VECTOR 
70  7  23 

WOULD    YOU    LIKE    TO    CHANGE    ANY     CF    THE    DATA? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 

ENTER    THE    NUMBER    OF    THE    YEAR     YOU    WISH    TO    SEE 
1  0 

DC    YOU     WISH     TC    SEE    THE    INTERVENING    YEARS? 

0  NO  - 

1  YES 


50 


T 

N 

PERCENT 

E 

0 

1 
2 
3 
TOTAL 

300  I 
200  < 
100  I 
600  I 

i        50) 

!     iooJ 

10 

1 

2 
3 

TOTAL 

351  | 

75  | 

302  | 

728  | 

48) 

10) 

41) 

121) 

100 


DO    YOU    WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

1 
ENTER    THE    NUMEER    OF    THE    YEAR    YOU    WISH    TO    SEE 

50 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
50  1  375    (  43) 

2  79     |  9) 


3  415     <  48 


IOTAL  868     (  145)  100 

DO    YOU    WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7  QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
100 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


100 


DC    YOU    WISH     TC    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7  QUIT    PROGRAM 

0 

DO    YOU     WISH     TC    SEE    THE    STEADY    STATE    VECTOR? 

0  NC 

1  YES 


100     1 

376  ( 

43 

2 

79  ( 

9 

3 

418  ( 

48 

TOTAL 

873  ( 

146 

1 
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999     1  376  (  43] 

2  79  (  9 

3  418  (  48 
TOTAL  873  (  100)        100 


AEE  YOU  THROUGH? 

0  NO 

1  YES 


1 
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APPENDIX  B 
COMPUTER  PRINTOUTS  OF  EXAMPLES  FOE  OPTION  2 

1 .   Exa  mple  (a) 

START 

DO  YOU  WISH  TO  ENTER  DATA? 

0  NC 

1  YES 

1 

ENTER  THE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    IENGTH    OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 

1 

ENTER    N  (INITIAL    STOCK    VECTOR) 

129    74    28    11 

ENTER    THE    PROMOTION    RATE    VECTOR. 

THIS    VECTOR     SHOULD    INCLUDE   THE    PROMOTION    RATES    FOR    THE    FIRST 
3    CLASSES. 

.  102    -046    .033 
ENTER    THE    WASTAGE    RATE    VECTOR. 

THE    4TH    VALUE    SHCULD    INCLUDE    THE    TOTAL    RATE    FOR       THAT 
STATE    DUE    TO    EITHER    PROMOTION    OR    WASTAGE. 

.  17    . 124    .1    .098 

ENTER  THE    NUMBER    OF    THE    RECRUIT    TYPE 

1  FIXED     RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 

• 

2 

ENTER    R     (RECRUITMENT    VECTOR) 

35000 

ENTER    ADDITIVE    INCREASE 


ENTER    THE   PERCENT    CODE 

0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7   QUIT  PROGRAM 

* 

1 
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WOULD  YOU  LIKE  TO  SEE  THE  ENTERED  DATA? 

0  NO 

1  YES 

1 

P    MATRIX 
0.728    0.102    0  0 

0  0.83       0.046    0 

0  0  0.867    0.033 

0  0  0  0-902 

N    VECTOR 
129    74    28    11 
OPTION    =2 
INC  =5 

TOTAI    RECRUITMENT 
35 

RECRUITMENT    PROPORTION    VECTOR 
10    0    0 

WOULD    YCU    LIKE    TO   CHANGE    ANY    OF   THE    DATA? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

« 

0 

ENTER    THE    NUMEER   OF    THE    YEAR    YOU    WISH    TO    SEE 
10 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 

T  N  PERCENT  R 


0            1 

129     < 

I                    53) 

2 

74     { 

31) 

3 

28     i 

12) 

4 

1  1     I 

5 

TOTAL 

242    | 

100) 

10            1  248  (  63) 

2  106  (  27) 

3  29  (                     7) 

4  10  (                      3) 

TOTAL  392  (  162)                     80 


DO    YOU     WISH     TO    SEE    ANY    OTHER     YEARS? 

0  NC 

1  Y  ES 

7  QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH     TO    SEE 
50 

DO    YOU     WISH     TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 


54 


50  1 

2 

3 

4 

TOTAL 


980 

523 

152 

39 

1695 


58 
31 
9 


2 

(  700) 


280 


DO    YOU    WISH    TO    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7         QOIT    PROGRAM 

1 

ENTER    THE    NUMEER   OF    THE    YEAR     YOO"    WISH    TO    SEE 
100 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 


530 


100     1 

1899  \ 

561 

31) 

2 

1075  i 

3 

343  | 

10) 

4 

102  { 

3 

TOTAL 

3419  { 

1413) 

DO    YOU    WISH    TO    SEE    ANY    OTHER     YEAPS? 

0  NC 

1  YES 

7         QUIT    PECGEAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SE3 
9  00 

DO    YOU    WISH    TO    SEE   THE    INTERVENING    YEAPS? 

0  NO 

1  YES 

0 


4530 


900     1 

16605  | 

54) 

2 

9898  | 

32 

3 

3395  I 

11 

4 

1130  | 

4) 

TOTAL 

31028  | 

12821) 

DC    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7  QUIT    PECGEAM 

* 

1 
ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
9  75 

DO    YOU     WISH     IC    SEE    THE    INTERVENING     YEARS? 

0  NO 

1  YES 

0 
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975            1  17984  <                    53] 

2  10725  (  32 

3  3681  (  IV 

4  1226  ( 
TOTAL  33617  (  13891)               4905 


DO    YOU     WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 

DO    YOU    WISH    TC    SEE    THE    STEADY    STATE    VECTOR? 

0  NO 

1  YES 

1 

PERCENTAGES     ARE    IN    STEADY    STATE 

S99  1  17984     (  53) 

2  10725  32) 

3  3681     (  11) 

4  1226     (  4) 

TOTAL  33617     (  100)  4905 

ARE    YOU    THROUGH? 

0  NO 

1  YES 


2.       Exa  mple     (b) 


START 

DO    YOU     WISH    TC    ENTER    DATA? 

0  NO 

1  YES 

1 

ENTER  THE    NUMEEP    Of    THE    MODEL    TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    LENGTH    OF    SERVICE 

3  MARKOV     GENERAL 

4  VACANCY 

3 
ENTER    N  (INITIAL    STOCK    VECTOR) 

300    200    100 

ENTER    P     (TRANSITION    MATRIX)      BY    ROWS 
ENTER     1  TH    ROW 

.6    .15     .05 

ENTER    2TH    ROW 

.75    .2    0 

ENTER    3TH    ROW 
.  05    0     .9 
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ENTER    THE    NUMEER   OF    THE    RECRUIT    TYPE 
1  FIXED    RICRUTT    VECTOR 


APpi±iyi_.(RECRUIT_SIZE) 
ADDITIVE       SYSTEM    SIZE) 


3  MULTIPLICATIVE     (RECRUIT    SIZE) 

— '  M 

5  MULTIPLICATIVE  '  (SYSTEM     SIZE) 

2 
ENTER    R     (RECRUITMENT    VECTOR) 

70   7    23 
ENTER    ADDITIVE    INCREASE 

25 

ENTER    THE    PERCENT    CODE 

0  NO    GRADE    PERCENTAGES 

1  GRADE    SIZE    AS    PERCENT   OF    TOTAL    GRADE    SIZE 

2  GRADE    SIZE    AS    PERCENT    OF    ORIGINAL    GRADE    SIZE 
7       QUIT    PROGRAM 

1 

WOULD    YCU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 

P    MATRIX 

0.6       0.15    0.05 
0.75    0.2      0 
0.05    0  0.9 

N  VECTOR 
300  200  100 
OPTION  =2 
INC     =25 

TOTAL  RECRUITMENT 
100 

RECRUITMENT  PROPORTION  VECTOR 
0.7  0.07  0. 23 

FOULD  YOU  LIKE  TO  CHANGE  ANY  OF  THE  DATA? 

0  NC 

1  YES 

7         QUIT    PRCGRAM 

0 
ENTER    THE    NUMEER   OF   THE    YEAR     YOU    WISH    TO    SEE 
10 
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DO    YOU    WISH    TC    SEE   THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 


T 

N 

PERCENT 

R 

0            1 
2 
3 

TOTAL 

300     | 
200     I 
100    { 

600     i 

50] 

33* 

17 

|                 100j 

1 

10            1 

2 
3 

TOTAL 

801     | 

162     | 

561     i 

1523     | 

53] 
IV 
37 
254; 

i 

325 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

1 
ENTER    THE    NUMEER   OF    THE    YEAR     YOU    WISH    TO    SEE 
50 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 
50  1 

2 
3 

TOTAL 


4374 
914 

4260 

9548 


1325 


DC    YOU    WISH     TC    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

1 
ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
100 

DO    YOU     WISH     TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


100 


0 


1 

2 
3 

TOTAL 


9066 

1903 

9469 

20438 


44] 

9' 

46] 

3406' 


DO    YOU     WISH     TC    SEE    ANY    OTHER     YEARS? 

0  NO 

1  YES 

7  QUIT    PROGRAM 


2575 


1 
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ENTER    THE    NUMBER    OF    THE    YEAR    YOU    WISH    TO    SEE 
3  50 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
350  1  32550    (  43) 

2  o853     {  9) 


3        35585  47 


TOTAL       74987  {     12498)      8825 


0 

400            1 

37246    ( 

43 

2 

7843     ( 

9 

3 

40808    ( 

48 

TOTAL 

85897     ( 

14316 

DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    FISH    TO    SEE 
4  00 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


10075 


DO    YOU     WISH    TC    SEE    ANY    OTHER    YEARS? 

0  N  0 

1  YES 

7  QUIT    PROGRAM 

0 

DO    YOU    WISH    TC    SEE    THE    STEADY    STATE    VECTOR? 

0  NC 

1  YES 

1 

PERCENTAGES     ARE    IN    STEADY    STATE 

999  1  37246     (  43] 

2  7843     (  9 

3  40808     j  48' 
TOTAL                 85897    (                  100)  10075 

ARE    YOU    THROUGH? 

0  NC 

1  YES 

1 
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APPENDIX  C 
COMPUTER  PRINTOUTS  OF  EXAMPLES  FOB  OPTION  3 

1  -   Example  (a ) 

START 

DO  YOU  WISH  TC  ENTER  DATA? 

0  NC 

1  YES 

* 

1 

ENTER  THE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

1 

ENTER    N  (INITIAL    STOCK    VECTOR) 

129    74    28    11 

ENTER    THE    PROMOTION    RATE    VECTOR. 

THIS    VECTOR     SHOULD    INCLUDE   THE    PROMOTION    RATES    FOR    THE    FIRST 
3    CLASSES. 

.  102    .046    .033 

ENTER    THE    WASTAGE    RATE    VECTOR. 

THE    4TH    VALUE    SHOULD    INCLUDE    THE    TOTAL    RATE    FOR       THAI 

STATE    DUE    TO     EITHER    PROMOTION    OR     WASTAGE. 

.17    .124    .1    .098 

ENTER  THE    NUMEER   OF    THE    RECRUIT    TYPE 

1  FIXED     RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 

3 
ENTER    R     (RECRUITMENT    VECTOR) 

35   0    0    0 
ENTER    MULTIPLICATIVE    FACTOR 

1  .01 

ENTER    THE    PERCENT    CODE 

0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7   QUIT  PROGRAM 

1 
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WOULD    YCa    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 

P    MATRIX 
0-728    0.102    0  0 

0  0.83       0.046    0 

0  0  0.867    0.033 

0  0  0  0.902 

N    VECTOR 
129    74    28    11 
OPTION    =3 
FACT  =1.01 

TOTAL    RECRUITMENT 
35 

RECRUITMENT    PROPORTION    VECTOR 
10    0    0 

WOULD    YOU    LIKE    TO    CHANGE    ANY    OF   THE    DATA? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 
ENTER    THE    NUMBER   OF    THE    YEAR     YOU    WISH    TO    SEE 
1  0 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
_T_  _N_  _PERCENT_      ■ R_ 

~0~      ~T~  "129  (~  ~53)~~ 

2  74  (  31( 

3  28  12] 

4  1  1  (                      5 
TOTAL  2  42  (  100) 


38 


DO    YOU    WISH     TO    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7  QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH     TO    SEE 
50 

DO    YOU     WISH     TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 


10     1 

137  | 

54) 
31 

2 

79  | 

3 

27  1 

11) 

4 

10 

V 

TOTAL 

253  i 

105 

61 


50     1 

204  < 

[       55) 
31 

2 

116  \ 

3 

37  | 

10) 

4 

11  I 

3 

TOTAL 

368  | 

152 

57 


DO    YOU    WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

• 

1 

ENTER    THE    NUMBER    OF    THE    YEAR    YOU    WISH    TO    SEE 
100 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 

100     1 

336  ( 

55 

2 

190  ( 

31 

3 

61 

10 

4 

19  ( 

6  0b 

3 

TOTAL 

250 

94 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

0 

DO    YOU    WISH    TO    SEE    THE    STEADY    STATE    VECTOR? 

0  NC 

1  Y  IS 

1 

PERCENTAGES    ARE    IN    STEADY    STATE 
999  1  336     (  55 

2  190     (  31 

3  61     (  10 

4  19    f  3 


TOTAL  606     {  100)  94 


o3| 


ARE    YOU    THROUGH? 

0  NO 

1  YES 

1 

2.       Exa  mple     (b) 

START 

DO    YOU     WISH     TO    ENTER    DATA? 

0  NC 

1  YES 
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ENTER  THE    NUMEER   OF   THE    MODEL    TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    LENGTH    OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 

3 

ENTER    N  (INITIAL    STOCK    VECTOR) 

300    200    100 
ENTER    P     (TRANSITION    MATRIX)     BY    ROWS 
ENTER    1  IH    ROW 

.6   .15    .05 

ENTER    2TH    ROW 

.75    .2    0 
ENTER    3TH    ROW 

.  05    0    .9 

ENTER  THE    NUMBER    OF    THE    RECRUIT    TYPE 

1  FIXED    RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 

3 

ENTER    R     (RECRUITMENT    VECTOR) 

70   7    23 

ENTER    MULTIPLICATIVE    FACTOR 

1  .01 

ENTER    THE    PERCENT    CODE 

0  NO    GRADE    PERCENTAGES 

1  GRADE    SIZE    AS    PERCENT    OF    TOTAL    GRADE    SIZE 

2  GRADE    SIZE    AS    PERCENT    OF    ORIGINAL    GRADE    SIZE 
7       QUIT    PROGRAM 

1 

WOULD    YOU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 

P    MATRIX 
0.6      0.15    0-05 
0.  75    0.2      0 
0.05    0  O.S 

N    VECTOR 
300    200    100 

OPTION    =3 

FACT         =1.0  1 

TOTAL    RECRUITMENT 
100 

RECRUITMENT    PROPORTION    VECTOR 
0.7    0.07   0-23 
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WOULD  YOU    LIKE    TO   CHANGE    ANY    OF    THE    DATA? 

0  NO 

1  YES 

7  QUIT    PROGRAM 

0 
ENTER    THE    NUMEER   OF    THE    YEAR    YOU    WISH    TO    SEE 
1  0 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 


PERCENT 


1 
2 

3 

TOTAL 


300 
200 
100 
600 


10 


1 
2 

3 

TOTAL 


370 

78 

312 

761 


33  i 
17 
100  i 

49 

10 

41 

127) 


109 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
50 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


50 


0 


1 

2 
3 

TOTAL 


575 
121 

603 
1298 


44) 
9 
46 
2  1  6  i 


163 


DO    YOU     WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

1 
ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
1  00 

DO    YOU     WISH    TC    SEE    THE    INTERVENING     YEARS? 

0  NO 

1  YES 

0 


64 


100 


1 

947 

2 

199 

3 

996 

TOTAL 

2141 

44) 
9 
47 
357' 


268 


DO    YOU    WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 

DO    YOU    WISH    TO    SEE    THE    STEADY    STATE    VECTOR? 

0  NO 

1  YES 

1 

PERCENTAGES  ARE    IN    STEADY    STATE 

999            1  947    <                    44) 

2  199    (  9) 

3  996     {  47) 

TOTAL  2141  100)  268 


ARE    YOU    THROUGH? 

0  NO 

1  YES 


1 
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APPENDIX    D 
COMPUTER   PRIHTOUTS   OF    EXAMPLES    FOR    OPTION    4 

1.      Case   J  Example    (a) 

START 

DO  YOU  WISH  TC  ENTER  DATA? 

0  NO 

1  YES 

1 

ENTER  THE  NUMBER  OF  THE  MODEL  TYPE 

1  MARKOV  HIERARCHICAL 

2  MARKOV  LENGTH  OF  SERVICE 

3  MARKOV  GENERAL 

4  VACANCY 

* 

1 

ENTER    N  (INITIAL    STOCK    VECTOR) 

129    74    28    11 
ENTER    THE    PROMOTION    RATE    VECTOR. 

THIS    VECTOR    SHOULD    INCLUDE   THE    PROMOTION    RATES    FOR    THE    FIRST 
3    CLASSES. 

.  102    .046    .033 

ENTER    THE    WASTAGE    RATE    VECTOR. 

THE    4TH    VALUE    SHOULD    INCLUDE    THE    TOTAL    RATE    FOR       THAT 

STATE    DUE    TO    EITHER    PROMOTION    OR    WASTAGE. 

. 17    . 124    . 1    .098 

ENTER  THE    NUMBER   OF    THE   RECRUIT    TYPE 

1  FIXED     RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 

4 
ENTER    RPROP      (RECRUITMENT    PROPORTION     VECTOR) 

10    0    0 

ENTER    ADDITIVE    INCREASE 

0 

ENTER    THE    PERCENT    CODE 

0  NO    GRADE    PERCENTAGES 

1  GRADE    SIZE    AS    PERCENT    OF    TOTAL    GRADE    SIZE 

2  3RADE    SIZE    AS    PERCENT    OF    ORIGINAL    GRADE    SIZE 
7       QUIT    PROGRAM 

1 

WOULD    YOU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 

€6 


P    MATRIX 
0.728    0.102    0  0 

0  0.83       0.  046    0 

0  0  0.867    0.033 

0  0  0  0.902 

N    VECTOR 
129    74    28    11 
OPTION    =4 
INC  =0 

RECRUITMENT    PROPORTION    VECTOR 
10    0    0 

WOULD    YOU    LIKE    TO   CHANGE    ANY    CF   THE    DATA? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
1  0 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 

T  N  PERCENT  R 


0            1 

129    | 

53) 

2 

74     | 

31) 

3 

28     l 

12) 

4 

11     < 

5 

TOTAL 

242    | 

100) 

10             1 

129     | 

I                    53l 
32) 

2 

77     { 

3 

27    | 

111 

4 

10     i 

4 

TOTAL 

242     1 

100) 

35 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7  QUIT    PRCGEAM 

1 
ENTER    THE    NUMEER   OF    THE    YEAR     YOU    WISH    TO    SEE 
50 

DC    YOU     WISH     TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  IES 

0 


50            1  129  (                    53) 

2  77  32) 

3  27  (  11) 

4  9  (  4) 

TOTAL  242  (  100)                     35 
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DO    YOU    WISH    TO    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

* 

1 

ENTER    THE    NUMBER    OF    THE    YEAR     YOU    WISH    TO    SEE 
100 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 


100            1 

129     | 

;              53) 

32 

2 

77     | 

3 

27    i 

11! 

4 

9     < 

4 

TOTAL 

242    ( 

100) 

35 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  Y£S 

7         QUIT    PRCGEAM 

0 

DO  YOU  WISH  TO  SEE  THE  STEADY  STATE  VECTOR? 

0  NC 

1  YES 


1 

PERCENTAGES     ARE    IN    STEADY    STATE 

999           1  129  (  53) 

2  77  j  32) 

3  27  J  11) 

4  9  (  4) 
TOTAL  242  (  100) 


35 


ARE    YOU    THROUGH? 

0  NC 

1  YES 


1 


Case     1 


Example    (b) 


START 

DC    YOU     WISH     TC    ENTER    DATA? 

0  NC 

1  YES 

1  • 

ENTER    THE    NUMEER    OF    THE    MODEL    TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    LENGTH   OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 
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ENTER    N  (INITIAL    STOCK    VECTOR) 

300    200    100 
ENTER    P     (TRANSITION    MATRIX)     EY    ROWS 
ENTER    1  1H    ROW 

.6.15     .05 

ENTER    2TH    ROW 

.75    .2    0 

ENTER    3TH    ROW 

.05    0    .9 

ENTER  THE    NUMEEE    OE    THE    RECRUIT    TYPE 

1  FIXED     RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 

4 
ENTER    RFROP     (RECRUITMENT    PROPORTION    VECTOR) 

.7   .07    -23 

ENTER    ADDITIVE    INCREASE 

0 

ENTER    THE    PERCENT    CODE 

0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7   QUIT  PROGRAM 

1 

WOULD    YOU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 


P    MATRIX 
0.6      0.15    0.  05 
0.75    0.2      0 
0.05    0  0.9 

N    VECTOR 
300    200    100 
OPTION    =4 
INC  =0 

RECRUITMENT    PROPORTION    VECTOR 
0.7    0.07   0.  23 


WOULD  YCU    LIKE    TO    CHANGE    ANY     OF    THE    DATA? 

0  NC 

1  YES 

7  QUIT    PROGRAM 

0 
ENTER    THE    NUMBER    OF    THE    YEAH     YOU    WISH    TO    SEE 
1  0 
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DO    YOU     WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 


PERCENT 


0            1 
2 
3 

TOTAL 

300     | 
200     { 
100     i 
600     \ 

50; 

33( 

17; 

;            io<r 

10            1 
2 
3 

TOTAL 

284     | 

61     i 

255     I 

600     I 

471 
10 

43 
100 

i 
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DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PHCGEAM 

1 
ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 

■ 

50 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 
50  1 

2 
3 


258 

54 

287 

600 


%3! 


48) 
TOTAL  600     (  100)  69 


DO    YOU    WISH     TC    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7  QUIT    PROGRAM 

1 
ENTER    THE    NUMEER    CF    THE    YEAR     YOU    WISH    TO    SEE 
1  GO 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 


100 


1 

2 
3 

TOTAL 


258 

54 

287 

600 


43] 

48 
100 


69 


70 


DO  YOU  WISH  TC  SEE  ANY  OTHER  YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 


DO    YOU    WISH    TO    SEE    THE    STEADY    STATE    VECTOR? 

0  NC 

1  YES 

m 

1 

PERCENTAGES    ARE    IN    STEADY    STATE 

999  1  258    (  43! 

2  54     (  9 

3  287     {  48 
TOTAL                      600     (                 100  69 


ARE    YOU    THROUGH? 

0  NC 

1  YES 

1 
3-       Case    2  Example    (a) 

START 

DO    YOU     WISH     TC    ENTER    DATA? 

0  NC 

1  YES 

1 

ENTER    THE    NUMBER   OF    THE    MODEL    TYPE 

1  MARKOV    HIERARCHICAL 

2  MARKOV    LENGTH    OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 

1 

ENTER    N  (INITIAL    STOCK    VECTOR) 

129    74    28    11 

ENTER    THE    PROMOTION    RATE    VECTOR. 

THIS    VECTOR     SHOUID    INCLUDE    THE    PROMOTION    RATES    FOR    THE     FIRST 
3    CLASSES. 

.  102    -046    .033 

ENTER    THE    WASTAGE    RATE    VECTOR. 

THE    4TH    VALUE    SHOULD    INCLUDE    THE    TOTAL    RATE    FOR       THAT 

STATE    DUE    TO     EITHER    PROMOTION    OR    WASTAGE. 

.  17    .  124    .  1    .098 

ENTER  THE    NUMEER    OF    THE    RECRUIT    TYPE 

1  FIXED     RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM     SIZE) 
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ENTER    RPROP     (RECRUITMENT    PROPORTION    VECTOR) 

10    0     0 

ENTER    ADDITIVE    INCREASE 


ENTER    THE    PERCENT    CODE 

0  NO    GRADE    PERCENTAGES 

1  GRADE    SIZE    AS    PERCENT   OF    TOTAL    GRADE    SIZE 

2  GRADE    SIZE    AS    PERCENT    OF    ORIGINAL    GRADE    SIZE 
7       QUIT    PROGRAM 

* 

1 

WOULD    YCU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

m 

1 

P    MATRIX 

0-728    0.102  0  0 

0  0-83  0.046    0 

0  0  C. 867    0.033 

0  0  0  0.902 

N    VECTOR 
129    74     28    1  1 

OPTION    =4 

INC  =5 

RECRUITMENT    PROPORTION    VECTOR 
10    0    0 

tfOULD    YOU    LIKE    IC    CHANGE    ANY     OF   THE    DATA? 

0  NO 

1  YES 

7         QUIT    PRCGRAM 

0 
ENTER    THE    NUMBER    OF    THE    YEAR     YOU    WISH    TO    SEE 
10 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 

T  N  PERCENT  R 


0     1 

129 

( 

53) 

2 

74 

} 

31) 

3 

28 

r 

12 

4 

11 

( 

5) 

TOTAL 

242 

( 

100)  _ 

10     1 

165 

( 

57} 
30) 

2 

89 

3 

28 

( 

10) 

4 

10 

I 

3) 

TOTAL 

292 

121) 

47 
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DO    YOa    WISH    TO    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PRCGRAM 

1 

ENTER    THE    NUMBER   OF    THE    YEAR     YOU    WISH    TO    SEE 

50 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 

50     1 

273  | 

56) 

2 

155 

31) 

3 

49  I 

10 

4 

15  ( 

3) 

TOTAL 

492  | 

203) 

76 


DO    YOU    WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7         QUIT    PRCGEAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAS     YOU    WISH    TO    SEE 
1  00 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 

100     1 

407  | 

55) 

2 

235  I 

32 

3 

77  | 

10) 

4 

24  i 

3 

TOTAL 

742  | 

307) 

113 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NO 

1  YES 

7  QUIT    PRCGEAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
9  00 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO      • 

1  YES 


0 

900     1 

2537  ( 

54) 

2 

1513  ( 

32) 

3 

519  ( 

1  1 

4 

173  ( 

4) 

TOTAL 

4742  ( 

1960) 

692 
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DO    YOU    HISH    TC    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
9  50 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 


950            1 

2670    | 

[                    53i 

2 

1593     i 

32) 

3 

547     i 

11 

4 

182     1 

4 

TOTAL 

4992    | 

2063) 

728 


DO    YCU    WISH    TO    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7  QUIT    PECG5AM 

0 

DO    YOU    WISH    TC    SEE    THE    STEADY    STATE    VECTOR? 

0  NC 

1  YES 


PEPCENTAGES     ABE  IN    STEADY    STATE 

999  1  2670     (  53' 

2  1593    (  32 

3  547     (  1V 

4  182     (  4j 
TOTAL                   4992     (                 100)  728 


ARE    YOU    THROUGH? 

0  NO 

1  YES 

1 

y-      Case   2  Example    (b) 

START 


DO    YOU     WISH     TC    ENTER    DATA? 

0  NC 

1  YES 

1 
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ENTER    THE    NUMEER   OF   THE    MODEL    TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    LENGTH   OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 


ENTER    N  (INITIAL    STOCK    VECTOR) 

300    200    100 
ENTER    P     (TRANSITION    MATRIX)     BY   ROWS 
ENTER    1  IH    ROW 

.6    .15    .0  5 

ENTER    2TH    ROW 

.75    .2    0 

ENTER    3TH    ROW 

.05    0    .9 

ENTER  THE    NUMEER   OF    THE   RECRUIT    TYPE 

1  FIXED    RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 


ENTER    RPROP     (RECRUITMENT    PROPORTION    VECTOR) 

.7   -07    .23 
ENTER    ADDITIVE    INCREASE 

25 

ENTER    THE   PERCENT    CODE 

0  NO    GRADE    PERCENTAGES 

1  GRADE    SIZE    AS    PERCENT   OF    TOTAL    GRADE    SIZE 

2  GRADE    SIZE    AS    PERCENT    OF    ORIGINAL    GRADE    SIZE 
7       QUIT    PROGRAM 

1 

WOULD    YOU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 

P    MATRIX 
0.6       0.15    0.05 
0.75    0.2      0 
0.05    0  0.9 

N  VECTOR 
300  200  100 

OPTION  =4 
INC     =25 

RECRUITMENT    PROPORTION     VECTOR 
0.7    0.07   0. 23 
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WOULD  YOU    LIKE    10    CHANGE    ANY     OF   THE    DATA? 

0  NO 

1  YES 

7  QUIT    PROGRAM 

0 
ENTER    THE    NUMEER    OF    THE    YEAR    YOU    WISH    TO    SEE 
1  0 

DO    YOU    WISH    TO    SEE   THE    INTERVENING    YEARS? 

0  NO 

1  YES 


0 


PERCENT 


1 

2 
3 

10TAL 


300 
200 
100 
600 


10 


1 
2 
3 

TOTAL 


417 

88 
345 
850 


50 

33 

17 

100 

49] 
10 
4V 
1421 


128 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7  QUIT    PROGRAM 


ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
50 


DO    YOU    WISH    IC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 


50 


1 
2 
3 

10TAL 


827 

173 

849 

1850 


45] 

9 

46 

303' 


239 


DO    YOU     WISH     TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

1 

ENTFR    THE    NUMBER    OF    THE    YEAR     YOU    WISH     TO    SEE 
1  00 

DO    YOU     WISH     TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
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100            1  1365  (  44) 

2  287  (  S) 

3  1448  {  47) 

TOTAL  3100  (  517                    382 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

1 

ENTER    THE    NUMEER   OF    THE    YEAR     YOU    WISH    TO    SEE 
3  50 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
350  1  4056     (  43) 

2  854     |  9) 


3  4440     (  47 


IOTAL  9350     (  1558)  1098 

DO    YOU    WISH    TC    SEE    ANY    OTHER    YEARS? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    HISH    TO    SEE 
4  00 


DO    YOU    WISH    TO    SEE   THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
400  1  4594     (  43) 

2  967     I  9) 


3  5039     (  48 


IOTAL  10o00     (  1767)  1241 


DO    YOU    WISH    TC    SEE    Ai*Y    OTHER    YEARS? 

0  NC 

1  YES 

7  QUIT     PRCGRAM 

0 

DO  YOU  WISH  TO  SEE  THE  STEADY  STATE  VECTOR? 

0  NC 

1  YES 

* 

1 
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PERCENT  AGES    AEE    IN    STEADY    STATE 
999  1  4594     ( 

2  967 

3  5039 
TOTAL                 10600 


43 
9 
48 
100) 


1241 


AEE    YOO    THROUGH? 

0  NO 

1  YES 


1 


78 


APPENDIX    E 
COMPUTER    PRIHTOOTS    OF    EXAMPLES    FOE    OPTION    5 

1.       Example     (a) 

S1ART 

DO    YOQ    WISH    TO    ENTER    DATA? 

0  NC 

1  YES 

1 

ENTER  THE  NUMEER  OF  THE  MODEL  TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    LENGTH   OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 

1 

ENTER    N  (INITIAL    STOCK    VECTOR) 

« 

129    7  4    28    11 

ENTER    THE    PROMOTION    RATE    VECTOR. 

THIS    VECTOR    SHOULD    INCLUDE   THE    PROMOTION    RATES    FOR    THE     FIRST 
3    CLASSES. 

.  102     .046    .033 

ENTER    THE    WASTAGE    FATE    VECTOR. 

THE    4TH    VALUE    SHOULD    INCLUDE    THE    TOTAL    RATE    FOR       THAT 

STATE    DUE   TO     EITHER    PROMOTION    OR    WASTAGE. 

.  17    .  124    . 1    .098 

ENTER  THE    NUMEER   OF    THE   RECRUIT    TYPE 

1  FIXED     RECRUIT    VECTOR 

2  AEDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 

4  ADDITIVE     (SYSTEM    SIZE) 

5  MULTIPLICATIVE     (SYSTEM    SIZE) 

5 
ENTER    RPROP     (RECRUITMENT    PROPORTION    VECTOR) 

10    0     0 
ENTER    MULTIPLICATIVE    FACTOR 

1  .01 

ENTER    THE    PERCENT    CODE 

0  NO  GRADE  PERCENTAGES 

1  GRADE  SIZE  AS  PERCENT  OF  TOTAL  GRADE  SIZE 

2  GRADE  SIZE  AS  PERCENT  OF  ORIGINAL  GRADE  SIZE 
7   QUIT  PROGRAM 

1 
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WOULD    YOU    LIKE    TO    SEE    TEE    ENTERED    DATA? 

0  NO 

1  YES 

1 

P    MATRIX 
0.728    0.102    0  0 

0  0.83       0.046    0 

0  0  0.867    0.033 

0  0  0  0.902 

N    VECTOR 
129    74    28    11 
OPTION    =5 
FACT  =1.01 

RECRUITMENT    PROPORTION    VECTOR 
10    0    0 

WOULD    YOU    LIKE    TC   CHANGE    ANY    OF   THE    DATA? 

0  NO 

1  YES 

7         QUIT    PROGRAM 

0 
ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
1  0 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 

_T_  _N_  _E!KCENT_  _R_ 

0           7  ^29  (  53] 

2  74  (  3V 

3  28  (  12] 

4  1  1  <  5 
TOTAL  242  (  100) 


10             1  148     (                     55 

2  82 

3  27 

4  10     I  4 


10     ( 
268     ( 


(  10) 


TOTAx.  268     (  111  42 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PRCGEAM 

1 

ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH     TO    SEE 
50 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 


30 


50     1 

222  | 

55) 
31 

2 

125  ( 

3 

40  { 

10) 

4 

12  | 

3) 

TOTAL 

400  I 

165 

62 


DO    YOU    WISH    TC    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PRCGfiAM 

1 

ENTER    THE    NUMEEfi   OF    THE    YEAR     YOU    WISH    TO    SEE 
1  00 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
100  1  362     (  55) 

2  206     (  31) 

3  66     (  10) 


4  20     (  3 


TOTAL  654     (  270  101 


DO    YOU    WISH    TO    SEE    ANY    OTHER    YEARS? 
0  NC 

1  V    P£ 

7  QUIT    PROGRAM 


DO    YOU    WISH  TO    SEE    THE    STEADY    STATE    VECTOR? 
C         NC 
1          YES 

1 

PERCENTAGES  ARE    IN    STEADY    STATE 

999           1  362     (                   55) 

2  2  06     (  31) 

3  06     {  10 

4  20     (  3) 

TOTAL  654     (                   100)                    10  1 


ARE    YOU    THROUGH? 

0  NO 

1  YES 

1 
2.       Example     (b) 

START 

DO    YOU     WISH     TC    ENTER    DATA? 

0  NC 

1  YES 

1 
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ENTER  THE    NUMEER   OF    THE    MODEL    TYPE 

1  MARKOV     HIERARCHICAL 

2  MARKOV    LENGTH    OF    SERVICE 

3  MARKOV    GENERAL 

4  VACANCY 

3 
ENTER    N  (INITIAL    STOCK    VECTOR) 

300    200    100 

ENTER    P     [TRANSITION    MATRIX)     BY    ROPS 
ENTER    1  TH    ROW 

.6.15    .OS 

ENTER    2IH    ROW 

.75    .2    0 

ENTER    3TH    ROW 

.  05    0     .9 

ENTER  THE    NUMEER   OF    THE    RECRUIT    TYPE 

1  FIXED     RECRUIT    VECTOR 

2  ADDITIVE     (RECRUIT    SIZE) 

3  MULTIPLICATIVE     (RECRUIT    SIZE) 
u  A  EDIT  IV  E     (SYSTEM    31  Zp) 

5  -IULTIPLICaTIVE     (SYSTE?     SIZE1 

* 

5 
ENTER  RFROP  (RECRUITMENT  PROPORTION  VECTOR) 

.7    .07    .23 

ENTER    MULTIPLICATIVE    FACTOR 

1  .01 

ENTER    THE    PERCENT    CODE 

0  NO    GRADE    PERCENTAGES 

1  GRADE    SIZE    AS    PERCENT    OF    TOTAL    GRADE    SIZE 

2  GRADE     SIZE    AS    PERCENT    OF    ORIGINAL    GRADE    SIZE 
7       QUIT    PROGRAM 

1 

WOULD    YOU    LIKE    TO    SEE    THE    ENTERED    DATA? 

0  NO 

1  YES 

1 

P    MATRIX 
0.6       0.15    0.05 

0.  75    0.2       0 
0.  05    0  0.9 

N    VECTOR 
300    200    100 
OPTION    =5 
FACT         =1.01 

RECRUITMENT    PROPORTION    VECTOR. 
0.7    0.07    0.23 

WOULD    YCU    LIKE    TO    CHANGE    ANY     OF    THE- DATA? 

0  NC 

1  YES 

7  QUIT    PROGRAM 

0 
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ENTER    THE    NUMBER    OF    THE  YEAR    YOU    WISH    TO    SEE 
1  0 

DO    YOU    WISH    TC    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 

_T  _N_                PERCENT_         _R_ 

0            T  300  (                 ~50) 

2  200  I  33) 

3  100  (  17) 

TOTAL  600  (                 100) 

10            1  31t>  (                    48] 

2  67  (  10 

3  277                        42 
TOTAL  66  1  {                 110                      87 


DO    YOU    WISH    TO    SEE    ANY    OTHER     YEARS? 

0  NC 

1  YES 

7         QUIT    PRCGEAM  " 

1 
ENTER    THE    NUMEER    OF    THE    YEAR     YOU    WISH    TO    SEE 
50 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

0 
50  1  435    (  4U) 

2  91     <  9) 


3  458 


TOTAL  985     (  1c4)  123 


46 
o4) 


DO    YOU     WISH     TC    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

1 

ENTER    THE    NUMEER    CF    THE    YEAR     YOU    WISH    TO    SEE 
1  00 

DO    YOU    WISH    TO    SEE    THE    INTERVENING    YEARS? 

0  NO 

1  YES 

• 

0 

100             1  717  (  44) 

2  150  (  9 

3  754  (  47] 
IOTAL  1621  (  270)                  20  3 
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DO    YOU     WISH    TO    SEE    ANY    OTHER    YEARS? 

0  NC 

1  YES 

7         QUIT    PROGRAM 

« 

0 

DO  YOU  WISH  TC  SEE  THE  STEADY  STATE  VECTOR? 

0  NC 

1  YES 

1 

PERCENTAGES     ARE  IN    STEADY  STATE 

999           1  717     (                   44] 

2  150    (  9 

3  754     (                   47 
TOTAL  1621     (                 100)                  203 

ARE    YOU    THROUGH? 

0  NC 

1  YES 

1 
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APPENDIX  F 
APL  COMPUTER  PBOGRAS  FOR  STEADY  STATE  DISTRIBUTION  FUNCTION 

This  is  the  APL  function  that  is  used  to  find  the  steady 
state  distribution  vector  for  all  options. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 


V    STDST ;  Q  ;  I  DENT ;  RV ;  W  ;  WSS  ;  NHAT1 ;  iVMr2  ;  F^CIT ;  IQ  ;  ZEP0  , 

Pt^(l,A')pPPP£P 

I  DENT*-  (K,K)p(l,KpO) 

W+(l,K)pl- (+/f 2 ]  P) 

G«-p+(  [>*/)  +  .  *Pf) 

+(TYPE=    1  2  3  u  5)/FIX,ADDREC,MULTREC/ADDSYS,MULTSYS 
FIX  :  A/*-P+  .  x  |  ( I  DENT-  P ) 

A'SS^/Vt(+/AO 

37-D  R0UND~N 

+END 
ADDREC : NHAT1+RPROP+ . *&{IDENT-P) 

NSS+NHAT1 * {+/NHAT1 ) 

+END 
MULTREC : FACTI+FACTxIDENT 

NHAT2+R+. *%{FACTI-P) 

NSS+NHAT2* ( +/NHAT2 ) 

+F.ND 
ADDSYS-.IQ+IDENT-Q 

IQ\ ;1]+Kpl 

ZERO+1  ,  (K-  DoO 

NSS+ZER0+. *a!Q 

-E77P 
MULT SYS : FACTI+FACTxIDENT 

NSS+(FACT-  1  )xRPR0P+.  ^{FACTI-Q) 
EWD1TCOUNT+999 

PERCENT+l (0 . 5+lOOxNSS) 

TOTPERCENT+100 

OUTPUT 
.    -»0 
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